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ON ELECTRON-CYCLOTRON WAVE RESONANCE IN 
INDUCTIVELY-COUPLED HF DISCHARGESt 


Edward P. Szuszczewicz* 

Laboratory for Planetary Atmospheres 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 

ABSTRACT 

The theory of Pft^'ffer on electron-cyclotron wave 
resonance phenomenorx , n inductively-coupled hf ring 
discharges has been extended by Sager to include electron 
temperature, non-uniform density distributions and the 
finite dimensions of a plasma of rectangular cross section. 
A scheme is herein proposed for adapting Sager's results 
to the experimentally important case of circular cross 
section. In order to test the scheme, an experiment was 
conducted in a low-pressure hf (27 MHz) argon discharge 
of cylindrical geometry with plasma densities in the range 
1.0 - 3.5(10 )/cm and static magnetic fields less than 
60 gauss. The agreement between theory and experiment is 
found to be good and lends credence to the scheme as a 
means for predicting resonance values of magnetic field 
in cylindrical hf discharges. 


t Part of this work was carried out while the author was 
supported by the Alexander von Humboldt Foundation as a 
guest scientist at the University of Wuerzburg (Germany) . 

* Presently a National Research Council Postdoctoral 
Research Associate. 


I . INTRODUCTION 


The superposition of a static magnetic field xoer- 

pendicular to the axis of a hf ring discharge n* ^ only 

optimizes carrier production but has also been observed to 

demonstrate a resonance-like behavior. These phenomena 

1 

were first reported by Neuert and theoretically explained 
2 3 

by Pfeiffer ’ as a result of a modified skin effect in an- 
isotropic plasmas and a resonant excitation of electron- 
cyclotron waves. 

2 4 

The theory of Pfeiffer has been extended by Sager 

to include a non-homogeneou*^ density distribution as well 
as temperature influences and the finite dimensions of a 
plasma of rectangular cross section. A scheme is herein 
proposed for adapting Sager's results to the experimentally 
important case of circular cross section and data are pre- 
sented which support the scheme. 

II. THEORETICAL CONSIDERATIONS 

To completely determine the electromagnetic field in 

a cylindrical hf discharge (hf magnetic field applied axially 

at the surface of the cylindrical plasma volume) with a 

static magnetic field superimposed in a direction normal 

to the cylinder axis is indeed very difficult. Pfeiffer's 

2 

approach to the problem was to restrict his considerations 
to a cold, homogeneous plasma slab which had its thickness 


2 


aligned with the superimposed static magnetic field. 

The electron-cyclotron resonance phenomenon was then theoreti- 
cally described by Pfeiffer by starting with the wave equation 
for the hf magnetic field H, Eq. (1), 

(7 X V X + €. . .) "S - 0 (1) 

o ij 

and requiring that the solution be linearly polarized in the 

plane of the plasma boundary. His solution to Eq. (1), where 

k is the vacuum wave number and e . , is the dielectric tensor 
o ij 

of the plasma, was a superposition of a standing right- and 
left-hand polarized wave which resulted from the Faraday 
effect. By imposing the condition that the hf electro- 
magnetic-field production of carriers in the plasma be at 
an optimum when the amplitude of the right-hand polarized 
wave (which propagates into the plasma along the static 
magnetic field B ) is large compared to that of the left-hand 
polarized wave, Pfeiffer obtained the relationship between 
the plasma density and the applied static magnetic field. 

For the ground resonance state this relationship is given 
in Eq. (2) 

•*res " ^®res " 

Where N is the plasma density, 3 is the electron-cyclotron 
frequency u)^ normalized to the angular frequency of the 


3 


generator uj, i.e. a « /'ju - (eB /m^)/<«, tii is the 

thickness oi the plasma slab and the subscript "res'* 

relers to the specific values under conditions of electron- 

cycloti’on resonance. The previously undefined quantities 

in Eq. (2) are the vacuum speed of light c, the permittivity 

of free space e^, and the mass and charge of an electron 

m and e, respectively, 
e 

4 

Like Pfeiffer, Sager began with the wave equation for 
the hf magnetic field (in a form different from that shown 
in Eq. (1)) but through a series of approximations has 
included the effects of non-zero plasma temperature, non- 
homogeneous density profiles, as well as a plasma geometry 
with rectangular cross-section. Since his work is relatively 
recent and has not been widely distributed a brief outline 
of his plasma model and theoretical considerations will be 
presented here. (It should be noted that the full scope of 
his approximation procedures and associated Justifications 
can best be appreciated by reading the original work.) 

Figure 1 presents the geometric configuration of Sager's 
plasma model where t|| is the plasma dimension parallel to 
the superimposed magnetic field and is the corresponding 
perpendicular dimension. (The dimension of the plasma 
normal to the plane of the figure is assumed to be infinite.) 
The net result of his approximation procedures is represented 
by Eq. (3) . 
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N 


rcb 


1 2 * 


(? -1*)^ t,i 

res II 
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(P 4 p^(\,)t*/— ) kT /m 

x’os i y t|| J e e 


(3) 

The temperature correction term appears in the denominator 

oi Eq. (3) and was introduced by Sager through the dependence 

of the plasma dielectric tensor on the electron temperature 

T . An approximate dependence of e . . on T was determined by 
e ij e 
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Sager by employing the results of the investigation of Mower 

on the conductivity of warm plasmas. 

The correction term for a non-homogeneous density distribu 

2 

tion appears in Eq. (3) as p (X^) and has been calculated by 
Sager in two limiting cases of ion mean-free-path X^. When 

X. is much greater than any plasma dimension t, Sager found 

2 2 
p (X. » t) =1.47; in the opposing limit he found p (x. « t) 

•la DL 

2.09. To arrive at these results Sager employed straight-line 
approximations (shown as dashed lines in Fig. 2) to the cosine- 
like profiles described in the early work of Tonks and Langmuir 
for the aforementioned limiting cases of ion mean-free-path. 
Shown as a solid curve in Fig. 2 is an experimental density 
profile observed in the cylindrical discharge which will be 
discussed in Section III. The experimental conditions were 
near resonance and X.^ was approximately equal to d/10, where 
d is the cylinder diameter. 

Of prime importance to the present investigation is the 
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Ion term which accounts lor the rectangular cross 
sv. ticiij () 1 the plasma volume. In Eq. (3) this term is 
designated by t* and is deiined in Eq. (4) in terms of the 
plasma cross sectional dimensions. 

t* (4) 

It was introduced by Sager in conjunction with his assumption 
Unit E<|. (5) represented the solution lor the right-hand 
polarized wave which propagated along the field B^. 

H - H Cos (k z) Cos (n -^) (5) 

ror r t 

JL 

In Eq. (5) k and H are the wave number and peak value, 
r or 

resijectively , of the right-hand polarized wave. The coordinate 
system is illustrated in Fig. 1 and defined so that y/\y\ 
and z/|J| II B^. 

To adapt the results of Sager's work to the experimentally 
important case of circular cross section, appropriate selectioni- 
must be made for t|| and t^ in terms of the cylinder diameter d. 

A reasonable approach is to have the circular cross section of 
a cylindrical discharge approximated by a square. With this 
assumption t|| = t^ and t = 1.41. 

The most important consideration revolves about the proper 
selection of tj| since in a circular cross section the thickness 
of the plasma in the direction parallel to B^ takes on continuous 
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IrcMi n iijj to a wuximuin whltii is oqual to the diameter 
d. This is illustrated In Fia* where it can readily be 
seen that t|()^ ■ - t,,)^ and that - t||)g * d. The 

scheme propgsed here, which di tiers from that oiJered by 

7 

Saj;er , involves a simple tiveraging of the values of ty over 
the circumference of the discharge. This averaging process 
is presented in Eq. (6) where the variable of integration 
is defined in Fig. 3 and where R is defined as the radius of 
the circular cross section. 


- / 

o 


y2 


2R Cos 0 d 


n 


( 6 ) 


The resulting adaption then approximates the circular 
cross section of diameter d by a square of side 2d/TT. 

It has been implicit in the scheme of adaptation that 
the excited hf waves ti’avei along the static magnetic field 
lines. This is certainly a reasonable consideration in the 
rectangular model of Sager for all waves initiated at the 
plasma boundary z = + t||/2. In the case of the circular 
cross section however the predominant wave-normal at the 
plasma boundary is aligned with the radius and forms an angle 
0 with . Now a characteristic of propagation in anisotropic 
plasmas is the fact that the direction of the wave normal 


is not the same direction in which the mean intensity of 
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*‘m‘r^y ilows. In u way analogous to the propagation ol 

B 

i'»iu,sphcTie whistlers the waves excited in the cylindrical 
i»l;isma volume will be "guided" by the magnetic tieid and 
will tend to travel along the lines of force. This will 
ho LL! listrated by employing the concept of a refractive 
index surface whose normal gives the direction of energy 
i 1 « »w . 


The conditions oi the experiments to be described in 

Section III are such that uj /ou > 1 where uj is the 

p res p 

electron plasma frequency. Neglecting electron-neutra I 
liollisions the Appleton-Hartree formula then gives Eq. (7) 
as the index of refraction for the single allowable mode of 


wave propagation. 


oj /u) 
P 


n - + 


Ii3 


cos 0-1 


(7) 


As be fore ft 


res 


^ ' res 

refers to the electron-cyclotron wave resonance 


value of r whose relationship to the plasma density is given 
Lu hq. (3). The refractive index surface corresponding to 
Eq. (3) is shown in Fig. 4 for B =2 and 3. Illustrated 
in the figure is an example of the energy flow direction with 
and without the plasma, s and s^ respectively, s and s^ are 
unit vectors with s being normal to the refractive index 

A 

surface and s^ being aligned with the radius vector. 


A A 


s*s^ = cos > and the positive sense of y is defined in the 
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C7(/^33*), shown as a dashed line, the waves can no longer 
penetrate within the plasma. These results point out that 
the physical situation is indeed quite conpatible with the 
adaptation scheme employed in the investigation and that the 
averaging process is justilied. 

III. EXPERIMENT 

In order to test the a lorementioned scheme ol adaptation 

a series ol measurements was carried out in an inductively- 

coupled hi discharge ol cylindrical geometry with length and 

diameter equal to 9.5 cm. The discharge was maintained by 

a single-turn coil (see Fig. 5) which was grounded at its 

point of symmetry and inductively coupled to a hi generator 

operating at 27 MHz. The plasma volume was bounded along its 

length by a metal plate on one end and by a wire grid on the 

other. The gas was argon and the operating pressure was 

-3 

maintained at 1.0 (10 ) torr. A pair of Helmholtz coils of 

mean diameter 46.8 cm provided the homogeneous static magnetic 
field which was varied in this experiment from 0 to approximately 
60 gauss. 


9 


The techuitiue employed in determining r , . was the 

3 

^ ;iiM* a^^ tha t (iriglnaily employed by Pieiiier and the re- 
sults are presented in Fig. 6 where iB the hi voltage 
iiieusured ;u*ro.ss the gup ol the single-turn secondaiy coil 

(Fiu. o) . N’ ^ has been previously delined and its value 
res 

was t*x per Imeu tally determined by employing the calculations 

9 

La 1 i‘umboise to analize the ion-saturation current collected 

b.v a Langmuir probe which was immersed in the plasma. The 

Langmuir probe was inserted along the axis ol the plasma 

<-y Under and was positioned in the geometrical center ol the 

plasma volume. The curves shown in Fig. 6 were generated by 

maintaining the middle-point plasma density constant and 

measuring the hi gap voltage as a function ol B. The minimum 

j.n the response ol U, . vs B then identifies B lor a 
^ hi res 

paiLlcular value ol N ^ . 

res 

Figure 7 includes two curves which result from Eq. (3) 

r 

as adapted in this paper to the circular cross section of a 

< vlindrical plasma volume. The curves have been calculated 

lor two diameters, 6.4 cm and 9.5 cm, under conditions which 

.set ^ . •'< d and T^ - 8(10‘^)®K. 
i 6 

The (N . ,, B results of Fig. 6 appear in Fig. 7 as 
solid dots and should be compared with the curve identified 
by d ^ 9.5 cm. The experimental conditions were such that 
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1 


✓J t> 

A mm i«j!u T ' 8V10 ) K and the eiTOi* burs represent 

1 e ^ 

the maximum possible inaccuracy in the experimental technique. 

In the experiments carried out in this investigation no 

value ol less than 2 was found lor values of N 

res res 

1 . 0 (10^^) /cm'^ . Oechsner^^ however measured resonance values 
Ox 8 less than 2 and his results are represented in Fig. 7 
by the open squai’es. His experiment was conducted in a 
cylindrical argon discharge ji diameter 9.5 cm with the 
pressure at 1.7(10“^) torr ard T - 6.5(10^) "k. It can 
be seen that the measurements of Oechsner and of this in- 
vestigation compare well with the appropriate theoretical 
curve (d == 9.5 cm) thus indicating that the scheme of adaptation 
yields a reasonable representation of the data for cylindrical 
hf dischrages. 

3 

As additional support of the scheme the data of Pfeiffer , 

-2 

which were collected in argon at 1.5(10 ) torr argon and 

which did not fit his original theoretical model, are plotted 
in Fig. 7 as open triangles and should be compared with the 
curve calculated for his experimental conditions at d = 6.4 cm. 

As in the case of the larger diameter, the agreement between 
the data and theory is good. 

IV COMMENTS & CONCLUSIONS 

A scheme has been proposed for adapting the theoretical 

4 

results of Sager to the experimentally important case of an 
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indue* tlvoly-<'ouplod cylindrical hf discharge. A prime con- 
sideration is the thickness of the plasma in the direction 
ol the superimposed static magnetic field. Since the field 
is applied in a direction perpendicular to the cylinder axis 
a simple averaging process over the circumference of the 
circular cross section has been employed. This process yields 
an average plasma thickness equal to 2d/rr and the net result 
of the proposed scheme approximates the circular cross section 
by a square of side 2d/rr. 

The experimental results of this investigation as well 

3 10 

as those of Pfeiffer and Oechsner agree well with the 
proposed scheme of adaptation and lend credence to the method 
as a means for establishing the relationship between plasma 
density and the static magnetic field in a cylindrical hf 
discharge under conditions of electron-cyclotron wave 
resonance . 
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H-IGURE CAPTIONS 


Figure 1 . 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Cross section of the hf plasma configuration 
assumed by Sager (Ref. 4) . 

Plasma density profiles along the superimposed 
static magnetic field. Dashed lines are those 
assumed in the plasma model of Sager (Ref. 4) 
and the solid curve is a profile observed in 
the cylindrical hf discharge studied in this 
investigation. 

An illustration of the variation of plasma 

— X 

thickness along the magnetic field which is 
superimposed perpendicularly to the axis of a 
cylindrical plasma volume. 

Refractive index surface. The length of the vector 

from the origin to the curve is nto/sju and 0 is 

P 

measured from the axis of the magnetic field . 
Schematic representation of the experimental arrange- 
ment employed in the study of the cylindrical hf 
discharge (shown in cross section). U, „ is the hf 

hi 

(27 MHz) voltage measured across the gap of the 
secondary coil. 
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Fif«ure 6. High-irequency voltage U. .. measured across 

ni 

the gap oi the secondary coil (Fig. 5) as a 

function of B, the electron-cyclotron frequency 

normalized to the generator frequency, at an argon 

-3 

pressure of 1.0(10 ) torr. The running parameter, 

N .^ ., is the middle-point plasma density which 

was maintained constant for the individual curves. 

Figure 7. Comparison between experimental results and the 

theory of Sager (Ref. 4) as adapted to the case 

of a cylindrical plasma volume. The data in this 

investigation have been taken from Fig. 6 and the 

theoretical curves have been calculated from 

Eq. (3) with p^(\.) “ 2.09, T * 8(10^.) °K, and 

1 e 

t\\ - - 2d/yi . The results of this investigation 

as well as those of Oechsner (Ref. 10) should be 
compared with the case d = 9.5 cm while the results 
of Pfeiffer (Ref. 3) should be compared with d = 6.4 cm. 
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